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BACKGROUND OF THE INVENTION 
Field of the Invention 

This invention relates to the use of Benders decomposition, a technique 
for solving multi-stage stochastic linear programming problems by decomposing them 
into a set of smaller linear programming problems, for supervisory control and 
coordination of independent tasks to be performed by a plurality of independent 
entities. More particularly, this invention relates to the application of Benders 
decomposition to the problem of scheduling maintenance in deregulated power 
systems. 

Description of Prior Art 

The electric power industry is moving toward a deregulated market- 
based operation in which the unbundling of services will establish independent entities 
such as power generators (GENCOs), transmission providers (TRANSCOs), 
independent customers (DISCOs) and independent system operators (ISOs). In the 
past, the first three entities were monopolized, owned and operated by regulated 
electric power companies. However, in the market-based structure, the supervisory 
control and coordination of these independent entities will be a vital task for power 
systems operators. 

In response to the electric industry restructuring, new phenomena, new 
circumstances, new risks, and new tools have emerged. Some of them arise due to the 
lack of experience with newborn issues, while others come as a necessity for the new 
structures. The electric power market differs from other commodity markets in that 
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electric power cannot be stored and must be consumed as generated. This condition 
drastically changes the nature of supervisory control and coordination in this industry 
as compared with that in other industries. As a result, applying theorems and models 
of other commodities to electricity have frequently misled participants of restructured 
electricity markets. In this regard, a systematic approach which will provide the 
proper mechanism as individual companies schedule their own activities, satisfy their 
own constraints and optimize their own objectives is desirable. From the customers 
perspective, maximizing the system availability and reliability and minimizing the 
cost of power delivery are desirable objectives which require coordination and control 
among these companies. 

Supervisory control may encompass different applications including the 
coordination of short-term and long-term maintenance scheduling of facilities in 
several independent power companies, fuel delivery and scheduling, and provision of 
ancillary services (dispersion of reserve capacity among power companies to sustain 
the continuity of service in the case of a contingency). It will be apparent that 
achievement of these objectives requires addressing various conflicting objectives, 
such as cost versus reliability, in power companies. It is imperative to maintain a 
certain level of reliability in a power system while individual companies within the 
system try to minimize their own cost of operation and maximize their revenues. 
Thus, it is clear that a systematic coordination and supervisory control is essential to 
assure the reliability of power delivery and to minimize power outages throughout the 
network. 
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Numerous methods relating to the planning and scheduling of 
maintenance for a variety of systems including power generation and distribution 
systems are disclosed by the prior art. For example, U.S. Patent 5,798,939 to Ochoa 
et al. teaches a computer workstation-based interactive tool for assessing the reliability 
5 of power systems, which tool can be used to determine the effect on the reliability of 
both substations and bulk generation and transmission systems of system additions, 
design alternatives, maintenance practices, substation configurations, and spare part 
policies. U.S. Patent 5,970,437 to Gorman et al. teaches a program for computerized 
management of plant maintenance which provides graphic representations of the 

10 mechanical and electrical systems of an operating plant where all of the components 
are treated as "objects" and associated through a relational database with text files 
providing component specifications, product identifications, component operating 
status, and all directly interconnected components including data as to flow direction 
and relative position in the system. For maintenance and service, each component is 

15 associated with existing and future specified work orders, and service requirements 

including consideration of probable life expectancy and the like. See also U.S. Patent 
5,311,562 to Palusamy et al. which teaches an integrated plant monitoring and 
diagnostic system for shared use by the operations, maintenance and engineering 
departments of a nuclear power plant; U.S. Patent 4,843,575 to Crane which teaches 

20 an interactive, dynamic, real-time management system comprising a plurality of 
powered systems and a central management facility; and U.S. Patent 6,006,171 to 
Vines et al. which teaches a computerized maintenance management system for the 
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process control environment which integrates a computerized maintenance 
management system with a process control system. 

SUMMARY OF THE INVENTION 
It is one object of this invention to provide a method for supervisory 
5 control and coordination of independent tasks in the deregulated electric power 

industry aimed at minimizing the possibility of blackouts, minimizing the cost of 
power delivery to customers, improving the system availability, and responding to 
environmental and regulatory concerns as individual power companies try to 
maximize their revenues. 

10 It is another obj ect of this invention to provide a method for supervisory 

control and coordination of independent tasks in the deregulated electric power 
industry which enables individual companies to schedule their own activities, satisfy 
their own constraints and optimize their own objectives while providing necessary 
coordination and control among these companies to maximize the system availability 

1 5 and reliability, and minimize the cost of power delivered to customers. 

These and other objects of this invention are addressed by a method for 
supervisory control and coordination of independent tasks to be performed by a 
plurality of independent entities comprising the steps of a) generating a plurality of 
schedules for performance of the independent tasks; b) submitting the plurality of 

20 schedules to a master coordinator for approval or disapproval of the schedules, 
resulting in generation of approval and/or disapproval decisions by the master 
coordinator; c) returning the approval and/or disapproval decisions to the independent 
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entities; d) adjusting the schedule for which a disapproval decision is returned, 
resulting in at least one adjusted schedule; and e) returning the adjusted schedule to 
the master coordinator for reconsideration. Steps b) through e) are repeated until all 
of the schedules have been approved. 

Benders decomposition with its mathematical feature is a particularly 
suitable match for providing assistance to human operators for supervisory control and 
coordination in the electric power industry in which there are a number of independent 
companies, generation, transmission, distribution, and customers, with their own 
objectives and there is a coordinator, for example the ISO in the power grid, which 
will coordinate the proposed activities. The application of Benders decomposition 
enables coordination of various conflicting objectives in different power companies, 
thereby enabling maintenance of a certain level of reliability in the power system 
while individual companies within the system try to minimize their own cost of 
operation and maximize their revenues. 

Although disclosed herein as being applicable to deregulated power 
systems, it will be apparent to those skilled in the art that the method of this invention 
employing Benders decomposition may be applied to any system requiring 
supervisory control and coordination of independent tasks to be performed by 
independent entities, and such applications are deemed to be within the scope of this 
invention. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

These and other objects and features of this invention will be better 
understood from the following detailed description taken in conjunction with the 
drawings wherein: 

Fig. 1 is diagram of an example of the Benders decomposition hierarchy 
in which the master problem represents the coordinator in the power grid and sub- 
problems identify individual power companies; 

Fig. 2 is a diagram of a Benders decomposition as applied to 
maintenance scheduling; 

Fig. 3 is a diagram showing application of the method of this invention 
as applied to multiple companies in a deregulated system; and 

Fig. 4 is a diagram of load data for an exemplary three-bus system. 

DESCRIPTION OF PREFERRED EMBODIMENTS 

Fig. 1 shows an example of the Benders decomposition hierarchy in 
which the master problem represents the coordinator in the power grid and sub- 
problems identify individual power companies. The individual companies schedule 
their own activities and submit their proposed schedules to the master problem which 
will check various criteria and approve or disapprove the proposed schedules. The 
coordinator's decision is returned as Benders cuts to sub-problems which will 
correspondingly adjust their schedules if their proposals are not approved. Subsequent 
iterations follow before the implementation of proposed schedules until all criteria are 
met. 
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In the past, when a power system was operated as a regulated monopoly, 
most techniques forthermal generator maintenance scheduling were based on heuristic 
approaches. These approaches consider a generating unit separately in selecting its 
optimal outage interval subject to constraints and an objective criterion for equalizing 
or leveling reserves throughout the planning interval, minimizing expected total 
production costs or leveling the risk of failure to meet demand. An example of a 
heuristic approach would be to schedule one unit at a time beginning with the largest 
and ending with the smallest. Most methods, mainly those based on heuristics, 
represent only the generation system and do not take into account the network 
constraint effects on the unit maintenance. The method of this invention extends 
Benders decomposition to include coordination between a GENCO and a TRANSCO 
for the inclusion of network constraints in the maintenance scheduling problem. 

The maintenance scheduling problem will determine the period for 
which generating units of a GENCO should be taken off line for plant preventive 
maintenance over the course of a one or two year planning horizon so that the total 
operating cost is minimized while system energy and reliability requirements as well 
as the number of other constraints in TRANSCO are satisfied. Since GENCO and 
TRANSCO are two separate entities, TRANSCO 's network constraints in 
maintenance scheduling of generating units may be excluded. The exclusion may 
result in an optimistic solution which will not satisfy network flow constraints and 
cannot be implemented since generating units in a GENCO are distributed throughout 
the network and interconnected by transmission lines. This may lead to different risk 
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levels to meet demand for a given amount of maintenance capacity outage depending 
on the unit location in the system. When network constraints are included, the 
problem becomes considerably more complex. The method of this invention is 
designed to address this problem. 

In coordinating the activities of GENCO and TRANSCO through an 
ISO, the objective is to minimize total operating and maintenance costs over the 
operating planning period, subject to unit maintenance and network constraints. To 
calculate the maintenance schedule, it is essential that numerous and complex 
constraints which limit the choice of scheduling times are incorporated into the 
solution method. Constraints in the maintenance scheduling problem are categorized 
as coupling and decoupling constraints. A coupling constraint is the requirement that 
generating units be overhauled regularly. This is necessary to maintain their 
efficiency at a reasonable level, maintain the incidence of forced outages low, and 
prolong the life of units. This procedure is incorporated periodically by specifying 
minimum/maximum times that a generating unit may run without maintenance. The 
time required for overhauling a unit is generally known, and hence the number of 
weeks that a unit is "down" is predetermined. It is assumed for purposes of this 
discussion that there is very little flexibility in manpower usage for maintenance. 
Furthermore, only a limited number of units may be serviced at one time due to the 
limited manpower. The available crew could be split into geographical and 
organizational types. For purposes of this description, we assume that the number of 
crew in each type required at each stage of overhaul of each unit is specified. 
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Network constraints in a TRANSCO in each time period are considered 



as decoupling constraints. The network can be modeled as either the transportation 
model or a linearized power flow model. For purposes of this description, the 
transportation model is used to represent TRANSCO's operating limits and the peak 
load balance equation. Mathematically, a complete generator maintenance schedule 
can be formulated as follows: 




S.T. 

maintenance constraints: 



x it =1 for / < e t or t > l t + d t 

x it = 0 for s t < t < s t + d t 

x it = 0 or 1 for e t < t < l t 



(i) 



1. crew availability 3. seasonal limitations 

2. resources availability 4. desirable schedule 
system constraints: 



(ii) 



(2-1) 



Sf+g+r = d Vf 
g< g.x Vt 
r<d Vt 



(hi) 
(iv) 
(v) 




(vi) 



(vii) 
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where: 



C u generation maintenance cost for unit i at time t 

c it generation cost of unit i at time / 

x lt unit maintenance status, 0 if unit is off-line for maintenance 

s t period in which maintenance of generating unit i starts 

e t earliest period for maintenance of generating unit i to begin 

l t latest period for maintenance of generating unit i to begin 

d t duration of maintenance for generating unit i 

r vector of dummy generators which corresponds to energy not served at time period t 

f maximum line flow capacity in matrix term 

/ active power flow in vector term 

g maximum generation capacity in vector term 

g vector of (gj power generation for each unit at time t 

d vector of the demand in every bus at time t 

S node-branch incidence matrix 

8 acceptable level of expected energy not served 

The unknown variable x it in (2-1) is restricted to integer values; on the 
other hand, g it has continuous values. Therefore (2- 1 ) corresponds to a mixed-integer 
programming problem. 

The objective of (2-1) is to minimize the total maintenance and 
production cost in a GENCO over the operational planning period. The first term of 
objective function (2-1) is the maintenance cost of generators; the second is the energy 
production cost. 

In (2-1), maintenance constraints (i) and (ii) represent GENCO's 
constraints. Also, constraints (i) represent the maintenance window stated in terms 
of the start of maintenance variables (s^. The unit must not be in maintenance before 
its earliest period of maintenance (e f ) and latest period of maintenance (e.g., l t +d^. 
The set of constraints (ii) consists of crew and resources availability, seasonal 
limitations, desirable schedule, and other constraints such as fuel and environmental 
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constraints. The seasonal limitations can be incorporated in and values of 
constraint (i). If units 1, 2 and 3 are to be maintained simultaneously, the set of 
constraints would be formed as follows: 

x u +x 2t +x 3t = 3 or x u +x 2t +x 3t = 0 

If in each maintenance area we have limited resources and crew 
available, the set of constraints would be formed as follows: 

Y.oJ\-x ll )<Z m (2-2) 

In the case of resource constraints, Z mt would be the amount of resource 
m available in area^4 for each time t and a mi would be a percentage of this resource 
required for unit i. In the case of crew constraints, the corresponding Z f would be the 
number of maintenance crew in area 0) and G mi would be a percentage of this crew 
required for maintenance of unit i. 

Constraints (iii)-(vi) represent the operation constraints checked by the 
ISO, in this case, peak load balance and other operation constraints such as generation 
and transmission capacity limits of the system. Constraint (vii) represents allowable 
energy unserved in the system. Benders decomposition is there applied to solve 
GENCO's and ISO set of equations in (2-1). 
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BENDERS DECOMPOSITION 

Before we discuss the solution methodology, let us present the Benders 

decomposition by considering the following general mixed integer program: 

Min Px+p(g) 
5 S.T. A l x<b l 

A^c+u{g)<b 2 (3-1) 
g>0 

Xf = 0 or 1 for all i 

where x is a vector of 0-1 variables with constant cost vector P and coefficient 

10 matrices A l and A 2 ; g is a vector of continuous variables with cost functions p and u; 

b l and b 2 are vectors of right hand side constants. Since the problem involves both 

discrete and continuous variables, it is unlikely that a direct approach to solve (3-1) 

would be computationally feasible. Instead the problem is partitioned as 

Min Px+Min g ^ 0 {p(g) | u(g) < b 2 - A^} 
15 S.T. A x x<b x 

x t = 0 or 1 (3-2) 
x € Q 

where Q is the set of x for which the constraints u(g)<b 2 -A 2 x can be satisfied. For 

each fixed x, the resulting inner minimization problem is 

20 Min p{g) 

S.T. u(g)<.b 2 -Ai>c (3-3) 

The Lagrangian relaxation of (3-3) is given by 

L(a ) = min s> _ 0 {p(g) + a (u(g) - (b 2 - A 2 x))} (3-4) 
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If g does satisfy w(g) < b 2 -A^c, the extra term in the objective will be non-positive and 
thus, for all Ct>0, 



The Lagrangian dual L is then defined by L = max a L(a). Under certain conditions 
sufficient for strong duality, 

L = min s> _ 0 { p(g) | u(g) < (b 2 - A 2 x)} ( 3-6 ) 

enabling us to replace the inner minimization of (3-2) by L. This replacement is 
justified later for our problem. With this replacement, (3-2) becomes 

Min Px + Max^ 0 {L ( O)} 
S.T. A lX <b, 

x t = 0 or 1 (3-7) 



L(a) < min s 



{piS) | KS) * (P 2 - A 2 x)} 



(3-5) 



! g>0 



If we let 



z =Px + max a ^ 0 {L( a) } 
= Px + max ai(i {mm gi0 { P (g) + a(u(g)-(b 1 -A 2 K))}} 



then (3-7) is equivalent to 



Min z 



S.T. A,x<b 



z > Px + min, 



{p(g)+ <Ku(g)-(b 2 -A 2 x))} for all a (3-8) 



Xf = 0 or 1 



xeQ 



Constraints (3-8) are referred to as feasibility cuts. 
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To complete the derivation of the master problem, the set of constraints that ensure 
xe Q must be characterized. This condition is satisfied if and only if 

max p> J\min % ^ {p(g) + fi(u(g) - (b 2 - A 2 x))}} < cc 

This condition is equivalent to 
5 min^ 0 {/3(u(g) - (b 2 - A 2 x))} < 0 for all J3 > 0 ( 3 " 9 ) 

Constraints of this form are referred to as infeasibility cuts. Thus our master problem 
is 

Minz 

S.T. z> Px+ min^ 0 {p(g) + a (u(g) - (b 2 - A 2 x))} for all a > 0 
: 10 mm^ 0 {p(u(g)-(b 2 - A 2 x))} <0 for all/? >0 (3-10) 

x t = 0 or 1 

The application of this method to our problem is discussed in the following section. 
SOLUTION METHODOLOGY 
Benders decomposition is applied to the generation maintenance 
15 problem as follows. If X denotes the vector of maintenance variables {x it }, Q 
represents the set of maintenance schedules for which constraints (iii)-(vii) are 
satisfied in all periods t, and operation cost w t is defined as 



then (2-1) can be written as 

Min X X C„ (1 - x„ ) + X Mm{w, \ (Hi) - {vii)} C 4 " 1 ) 
S.T. 



20 
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maintenance constraints: 

x it = 1 for t < e i or t> l t + d t 

x it = 0 for s t < t <s t + d t (i) 

5 x it = 0 or 1 for e t < t < l t 

1. crew availability 3. seasonal limitations 

2. resources availability 4. desirable schedule (ii) 

X€Q 

If the t th sub-problem was a linear program, it could be replaced by its dual as is done 
10 in the standard Benders decomposition. The Lagrangian dual of the f h sub-problem 
is given by 

where L t (K,H-,Y,C,|-0 is the Lagrangian function and K,7I, Y>CM- and [I are multipliers 
of constraints (iii)-(vii). 



15 



20 



25 



(4-3) 

The f th sub-problem is then placed by L t , and (4-1) is rewritten as 

M«IIC,(l-,J + 4 (4-4) 

t i 

S.T. 

maintenance constraints: 

x u = 1 for t < e; or t> l { + d i 

x it = 0 for s t < t < s t + d t (i) 

jc ; - ( = 0 or 1 for e i < t < l t 
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1. crew availability 3. seasonal limitations 

2. resources availability 4. preschedule (ii) 

XE Q 

To ensure X E Q, the maintenance schedule must ensure that sufficient reserve exists 
to provide a secure supply while minimizing the cost of operation. The f 1 sub- 
problem is feasible if and only if the optimal value of the following problem is less 
than e 

Min ^ r u 

S.T. Sf + g + r = d (4-5) 
r<d 

l/N/ 

Its Lagrangian dual is 

max vA ^ 0 U t {v,X,r,Tj) 

where U t (v,A,T,T|) is the following Lagrangian function and v,\,l and T] are 
multipliers of constraints (iii)-(vii). 



U,{v,Z,t,t]) = min g > 0 < 



I r* + 1 v„ ((X S lk f h ) + g, ( +r It -d,?j +1 X u (g lt -g .*.) + 



We then arrive at the generalized Benders master problem: 

Min z 
S.T. 

z>^C, ( (l-^) + Il,(w^) for all Ar,*,r,£//*0 

X U t (y,A,T,7j)<e for all v,A,r,r/>0 
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maintenance constraints: 




(4-6) 



1. crew availability 3. seasonal limitations 

2. resources availability 4. desirable schedule 



(ii) 



The problem is decomposed into a master problem and a set of independent operation 
sub-problems. The master problem, which in this model is an integer programming 
problem, is solved to generate a trial solution for maintenance schedule decision 
variables. This master problem is a relaxation of the original problem in that it 
contains only a subset of constraints. Its optimal objective value is a lower bound on 
the optimal value of the original problem. Once x lt variables are fixed, the resulting 
operation sub-problem can be treated as a set of independent sub-problems, one for 
each time period t, since there are no constraints across time periods. The set of 
operation sub-problems are then solved using the fixed maintenance schedule obtained 
from the solution of the master problem. At each iteration, the solution of sub- 
problems generates dual multipliers which measure the change in either production 
cost or reliability resulting from marginal changes in the trial maintenance scheduling. 
These dual multipliers are used to form one or more constraints (known as cuts) which 
are added to the master problem for the next iteration. The problem continues until a 
feasible solution is found whose cost is sufficiently close to lower bound, as shown in 
Fig. 2. 
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The initial maintenance master problem is formulated as follows: 
Min z 

S.T. z>XI{C„ (1-xj) 

maintenance constraints: 

x it =1 for t < e t or t > /,. + d t 

x it = 0 for s t < t < s t + d t 

x it = 0 or 1 for e t < t < l t (4-7) 

1. crew availability 3. seasonal limitations 

2. resources availability 4. desirable schedule 

Operation Sub-problem 

The sub-problem may not have any solutions due to the fact that the 
unserved energy cannot be kept above a desired level. If a sub-problem is infeasible, 
then an infeasibility cut is generated. For each infeasible sub-problem resulting from 
the n th trial solution of the master problem, the infeasible cut is of the form 

S^+S^ft^-^J^ (4-8) 

The multiplier A" may be interpreted as a marginal decrease in energy not supplied 

with a 1 MW increase in generation, given the n th trial maintenance schedule. The 
infeasibility cuts (4-8) will eliminate maintenance values which are not possible to be 
scheduled. 

If the sub-problem is feasible, then the fuel cost for period t, w t , depends 
on the utilization of the available units to satisfy load constraints in each time period 
subject to maintaining the unserved energy above a certain level. Thus the generation 
cost in period t can be expressed as 
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w t = Min c *git 
S.T. Sf+g + r = d 
g^g x" 

< 

The solution of the sub-problem is not complicated, since knowing 
which generators and transmissions are available during period / allows us to minimize 
10 the operation cost. The feasible cut is of the form 

z > x [ w ; + X (q a -xn)+ < & (*: - x n ;)) (4_ 10) 

where w" is the expected fuel cost for period / associated with the n th trial solution. 

The multiplier ?r" may be interpreted as the marginal cost associated with a 1 MW 

1 5 decrease in the power capacity, given the n th trial maintenance schedule. The cost cuts 
(4-10) will tend to increase the lower bounds obtained from each successive 
maintenance sub-problem solution. 

Maintenance Master Problem 
The maintenance master problem minimizes maintenance cost subject 

20 to maintenance constraints as well as feasibility and infeasibility cuts from the 
operation sub-problems. If all sub-problems are feasible, then their solutions yield a 
set of dual multipliers from which a feasibility cut is constructed. If one or more 
operation sub-problems are infeasible, then, for each infeasible sub-problem, an 
infeasibility cut is generated. 
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Min z 
S.T. 

z>SS{C„ (1-xj) 

t I 

• maintenance constraints: 

x it = 1 for t < e i or t > /,- + d t 

x it = 0 for s t < t < Si + df 

5 x it = 0 or 1 for e t < t < l t 

1. crew availability 3. seasonal limitations 

2. resources availability 4. desirable schedule 

• feasibility and infeasibility cuts from previous iterations 
10 • if all sub-problems are feasible then the feasible cut is: 

2> z (< + s c„ a -*"u)+ w - *u )) 

• if one or more sub-problems are infeasible, then the infeasible cuts are: 

15 

2 r" t + * g } (x" t - x u ) < s V? e infeasible sub-problem 

x€0,l (4-11) 

where 

n is the current number of iterations 

20 sfr, 7? are the multiplier vectors at n* iteration 

The important feature of the Benders decomposition is the availability 
of upper and lower bounds to the optimal solution at each iteration. These bounds can 
be used as an effective convergence criterion. The critical point in the decomposition 
25 is the modification of objective function based on the solution of the operation sub- 
problem. Associated with the solution of the operation sub-problem is a set of dual 
multipliers which measure changes in system operating costs caused by marginal 
changes in the trial maintenance. These multipliers are used to form a linear 
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constraint, written in terms of maintenance variable x. This constraint, known as 
Benders cut, is returned to the maintenance problem which is modified and solved 
again to determine a new trial maintenance plan. 

EXAMPLE 

We use a three-bus system as an example. The maximum energy not 
served requirement (e) is 0 p.u., and generator and line input data in per unit for 
GENCO and TRANSCO are given in Tables 1 and 2. Load data are depicted in Fig. 2. 
GENCO will perform maintenance on at least one generator. We assume the study 
period represents one time interval. Loads are assumed constant during the study 
period. 



Table 1. Generator Data (GENCO) 





Min Capacity 


Max Capacity 




Maint. Cost/Unit 


Unit 


(p.u.) 


(p.u.) 


Cost ($) 


($) 


1 1 


0.5 


2.5 


10 g! 


300 


2 


0.6 


2.5 


10 g 2 


200 


1 3 


0.6 


3.0 


10 g, 


100 



Table 2. Line Data (TRANSCO) 

Line Q/line # of lines Capacity/line (p.u.) 



1- 2 0.2 2 0.25 

2- 3 0.25 2 0.5 
1-3 0.4 2 025 



First, we solve the initial maintenance master problem. 
GENCO 's maintenance problem (Iteration 1): 
Min z 
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S.T. 300*(l-x 1 )+200*(l-x 2 )+100*(l-x 3 ) < z 
x l +x 2 +x- 3 < 2 
Xj < 1 x 2 < 1 x 3 < 1 

5 The solution is: x { =l x 2 =l x 3 =Q z=100 

ISO's operation problem (iteration 1): 

We check the feasibility of the operation sub-problem given the first trial 

of maintenance schedule. The feasibility check is given as follows: 

Min r x +r 2 +r 3 

10 S.T. -/i2-/i3 + gi +r i = l Load balance at bus 1 

-^3 + ii2 + ^2 +r 2 =3 Load balance at bus 2 

f u +f 2 3 + g2 + r 3 =\ Load balance at bus 3 

0.5<g!<2.5 Generator 1 limit 

0.6<g 2 <2.5 Generator 2 limit 

15 0.0<^ 3 <0.0 Generator 3 limit 

-2*0.25<f 12 <2*0.25 Line 1-2 flow limit 

-2*0.25 <f 13 <2*0.25 Line 1-3 flow limit 

-2*0.5 <f 23 < 2*0.5 Line 2-3 flow limit 

The primal solution of the feasibility check is: 

20 r=0.5 gl =2 g 2 =2.5 g,=0 f y2 =0.5 / 13 =0.5 / 23 =0 

The dual price of the operation sub-problem is: A, gl =0 ^- g 2 = l ^ g 3 = l 
The above LP solution is infeasible, since r 1 +r 2 +r 3 >0. The generation cost is set 
arbitrarily to 1000 because the solution is infeasible. The infeasible cut is: 
0 .5+ 1 *2. 5( 1 -x 2 )+ 1 *3 *(0-x 3 ) < 0 

25 GENCO's maintenance problem (iteration 2): 
Min z 

S.T. 300*(l-x 1 )+200*(l-x 2 )+100*(l-x 3 ) < z 

0.5 + 1*2.5* (1-jc 2 ) + 1*3* (0 - x 3 ) < 0.5 

x ] +x 2 +x i < 2 
30 Xj<1 x 2 <l x 3 <l 
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The solution is: x x =l x 2 =0 x 3 =l 



z=200 



ISO's operation problem (iteration 2): 



The feasibility check is as follows: 



Min 



S.T. 



-f23 + fn + g2 +r 2=3 
fl3 + f23 + g3 +r 3 =1 



0.5<g!<2.5 
0.0<g 2 <0.0 
0.6<g 3 <3.0 



Load balance at bus 
Load balance at bus 
Load balance at bus 
Generator 1 limit 



1 
2 
3 



Generator 2 limit 
Generator 3 limit 



-2*0.25 <f 12 <2*0.25 
-2*0.25<f 13 <2*0.25 
-2*0.5 <f„< 2*0.5 



Line 1-2 flow limit 
Line 1-3 flow limit 
Line 2-3 flow limit 



The primal solution of feasibility check is: r= 1.5 g t =1.5 g 2 =0 g 3 =2 / 12 -0.5 



The above LP solution is infeasible, since r t +r 2 +r 3 >0. The generation cost is set 
arbitrarily to 1000 because the solution is infeasible. The infeasible cut is as follows: 
1.5+l*2.5*(l-x 2 )<0 

GENCO's maintenance problem (iteration 3): 
Min z 

S.T. 300*(l-x 1 )+200*(l-xg 2 )+100*(l-xg 3 ) < z 
0.5 + 1*2.5* (l-x 2 ) + 1*3* (0 - jc 3 ) < 0.5 
1.5 + 1*2.5* (l-x 2 ) < 0 

x l +x 2 +x 3 < 2 

X,<1 JC 2 < 1 X 2 <1 

The solution is: x { =0 x 2 =l x 3 =l z=300 



/i3=0 / 23 -l 



The dual price of the operation sub-problem is: ^ g i=0 ^ g2 = l ^ g 3 - 0 
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Given the trial maintenance schedule in iteration 3, we apply the feasibility check as 



10 



15 



20 



25 



before which gives r=0. This means the trial schedule is feasible now. 
ISO's operation problem (iteration 3): 
The feasible problem is as follows 



Min w=10*g 2 +10*g 3 +300 

S.T. -fvfu + gi + r^l 

-f23 + fn+g2 + r 2 = 3 

/l3+/ 2 3+£ 3 +r3=l 

O.O^g^O.O 
0.6<g 2 <2.5 
0.6<g 3 <3.0 



-2*0.5 <f„ < 2*0.5 



Load balance at bus 
Load balance at bus 
Load balance at bus 
Generator 1 limit 
Generator 2 limit 
Generator 3 limit 
Line 1-2 flow limit 
Line 1-3 flow limit 
Line 2-3 flow limit 



The primal solution is: w=350 g 2 = 2-5 #3=2. 5 / 12 =-0.5 fn = -0.5 f 2 j,~^ 
The dual price of the operation sub-problem is: 7l gl =0 7l g2 =-10 7t g 3=-10 
The feasible cut for the third iteration is: z>350-10*2.5*(l-x 2 )-10*2.5*(l-x 3 ) 
GENCO's maintenance problem (iteration 4): 
Min z 

S.T. 300*(l-jc 1 )+200*(l-x 2 )+100*(l-x 3 ) < z 

0.5 + 1*2.5* (l-x 2 ) + 1*3* (0 - x 3 ) < 0.5 
1.5 + 1*2.5* (1-jc 2 ) < 0 
z>350-10*2.5*(l-x 2 )-10*2.5*(l-x 3 ) 

x ] +x 2 +x- i < 2 
x x <\ x 2 <\ x 3 <l 

olution is: Xj=0 x 2 

We stop here since z=w which means the cost is equal to the lower bound. 
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While in the foregoing specification this invention has been described 
in relation to certain preferred embodiments thereof, and many details have been set 
forth for purpose of illustration, it will be apparent to those skilled in the art that the 
invention is susceptible to additional embodiments and that certain of the details 
5 described herein can be varied considerably without departing from the basic 
principles of the invention. 
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